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We study the acceleration of heavy nuclei at SNR shocks when the process of ionization is taken 
into account. Heavy atoms {Zn > few) in the interstellar medium which start the diffusive shock 
acceleration (DSA) are never fully ionized at the moment of injection. The ionization occurs 
during the acceleration process, when atoms already move relativistically. For typical environ- 
ment around SNRs the photo-ionization due to the background galactic radiation dominates over 
Coulomb collisions. The main consequence of ionization is the reduction of the maximum energy 
which ions can achieve with respect to the standard result of the DSA. In fact the photo-ionization 
has a timescale comparable to the beginning of the Sedov-Taylor phase, hence the maximum en- 
ergy is no more proportional to the nuclear charge, as predicted by standard DSA, but rather to 
the effective ions' charge during the acceleration process, which is smaller than the total nuclear 
charge Z^- This result can have a direct consequence in the prediction of the knee structure of the 
cosmic ray spectrum. Moreover the acceleration of ultra-heavy elements beyond the Iron's maxi- 
mum energy is very hard to achieve making unlikely their possible contribution to the cosmic ray 
spectrum in the transition region between Galactic and extragalactic component. 
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1. Introduction 



The bulk of Galactic cosmic rays (CRs) is largely thought to be accelerated at the shock waves 
associated with supernova remnants (SNRs) through the mechanism of diffusive shock accelera- 
tion (DSA). A key aspect of DSA is the maximum energy achieved by different nuclei which is 
intimately connected with the interpretation of the knee structure and of the transition region from 
Galactic to extragalactic CR component. The knee is commonly interpreted as due to the super- 
position of the spectra of all chemicals with different cutoff energies. Using the flux of different 
components measured at low energies it has been shown that the knee is well reproduced if one 
assumes that the maximum energy of each specie £'max,A' is proportional to the nuclear charge Z^r 
[11]. Nevertheless the superposition of subsequent cutoff seems to be confirmed by the measure- 
ments of the spectrum of single components in the knee region: data presented by the KASCADE 
experiments show that the maximum energy of He is ~ 2 times larger than that of the protons [^. 
The relation £'max,A' °^ '^n is clearly predicted by DSA if one assume that the diffusion coefficient 
is rigidity-dependent and that nuclei are completely ionized during the whole acceleration process. 
Indeed the second assumption does not hold in general. When ions are injected into the accelera- 
tion process they are unlikely to be fully ionized, especially if of high nuclear charge. In fact the 
temperature of the circumstellar medium where the forward shock propagates varies from 10^ up to 
10^ K. Ifr^lO^K even hydrogen is not fully ionized, as demonstrated by the presence of B aimer 
lines associated with shocks in some young SNRs [^. For T ~ 10^ only atoms up to Zat = 5 can 
be completely ionized. The ordinary assumption made in the literature is that atoms are completely 
stripped soon after the beginning of the acceleration process. In spite of this assumption in [12] we 
showed that, for a typical SNR shock, the ionization time is comparable with the acceleration time, 
hence electrons are stripped long time after the injection, when ions already move relativistically. 
This fact has two important consequences: 1) the maximum energy of ions can be reduced with 
respect to the standard prediction of DSA and 2) the ejected electrons can easily start the acceler- 
ation process providing a source for the synchrotron radiation. The latter point has been already 



analyzed in [ ]12[ ] and [|13[]. Here we want to describe how the ionization can affect the maximum 
energy achieved by heavy ions during the acceleration at SNR shocks. 



2. Ionization vs. acceleration time 



In order to get the maximum energy achieved by different chemical specie, we need to com- 
pare the acceleration time with the ionization time. Let us consider atoms of a single specie N 
with nuclear charge Zf^ and mass m^r = Anip, which start the DSA with initial charge Z <Zj^ and 
momentum pinj. We adopt the acceleration time as computed in the framework of linear DSA for 
plane shock geometry: 

f'^(2M^'m) ^JL^,,,,»Jizlsls-w {^) yr. (2.1) 

Jp,n,Ul-U2\ Ml U2 J p niNC ° \Z) 

Here u is the plasma speed in the shock rest frame, and the subscript 1 (2) refers to the upstream 
(downstream) quantities (note that Wshock = u\). The downstream plasma speed is related to the 
upstream one through the compression factor r, i.e. uj =u\lr. We limit our considerations to 
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Strong shocks, which have r = 4. The turbulent magnetic field responsible for the particle diffusion 
is assumed to be compressed downstream according to B2 = rB\. The last equality in Eq. ( ^j] ) 
holds because we assume Bohm diffusion coefficient, that is Db = ripc/S, where jSc is the particle 
speed and = pc/ZeB is the Larmor radius. Here the magnetic field is expressed in /iG and the 
shock speed is ui = wslO^^cm/s. 



Eq. (2. 1 ) has to be compared with the ionization time. Ionization can occur either via Coulomb 



collisions with thermal particles or via photo-ionization with background photons. In |13] we 
showed that Coulomb collisions are negligible when the density of the circumstellar medium is 
«i < 30cm^^. This condition is satisfied for almost all SNR environments, with the exception of 
the very initial phase of the explosion, when the shock propagates in the dense progenitor's wind, 
or when the remnant collides with a molecular cloud. Here such scenarios will be neglected. 

The full ionization time results from the convolution of the photo-ionization cross section with 
the photon energy spectrum, driph/de, due to the total interstellar radiation field (ISRF), i.e.: 

^pi!ir) = J de^^^ca,i,{Ye), (2.2) 

where e' = ye is the photon energy as seen in the ion rest frame, which moves with Lorentz 
factor 7. The photo-ionization cross section can be estimated using the simplest approxima- 
tion for the ^-shell cross section of hydrogen-like atoms with effective nuclear charge Z, i.e. 
0'ph(£') = 64 a^^ OtZ^^ (In^z/s'^^^ where Ot is the Thompson cross section, a is the fine 
structure constant and I^^z is the ionization energy threshold for the ground state of the chemical 
specie N with Z^ — Z electrons. The numerical values of In^z can be found in the literature (see e.g. 



3. Maximum energy of Ions 

We stress that a consistent treatment of the ionization effects would require the use of time- 
dependent calculation. However, in order to get an approximate result for the maximum energy we 
can use the quasi-stationary version of the linear acceleration theory. 

First of all we assume that the maximum energy is achieved at the beginning of the Sedov- 
Taylor phase (tsr) In fact, for later times, t > tsr, the shock speed decreases faster than the 
diffusion velocity, hence particles at the maximum energy can escape from the accelerator and the 
maximum energy cannot increase further. If the total ionization time is comparable, or even larger 
than tsT, we do expect that ^max < ^^max' where we call E^^^ the maximum energy achieved by 
ions which are completely ionized since the beginning of acceleration. We also assume that initial 
level of atoms' ionization is only determined by the circumstellar medium temperature. We neglect 
complications arising from the dust sputtering process which could play an important role in the 
injection of ions in the DSA mechanism [Q]. 

Now let us consider a single ion injected with momentum pinj and total charge Zi. The ion 
undergoes acceleration at a constant rate oc Z\ during a time equal to the ionization time needed to 
lose one electron, Tphj, when it achieves the momentum pi. Tph.i and p\ can be determined simul- 
taneously equating the acceleration time with the ionization time, i.e. Tph,i(pi/niArc) = Tacc{po,Pi)- 
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Using Eq. (2. 1 ) the last condition gives: 

P\ _ Pin] ZiBiu\ 



Tph,i (pi/'wyvc) , (3.1) 



nipC nipC 1 .7 Zf^ 

where the subscript 1 label the quantities during the time Tphj. Because the photo-ionization occurs 



when ions akeady move relativistically [[l2|], we set j3 = 1. In Eq. ( |37T| ) B is expressed in jxG and 



Ush in units of 10 cm/s, while Tph.i is expressed in yr. Once the background photon distribution is 



known Eq. ( [3.1| ) can be solved numerically. After the first ionization the acceleration proceeds at a 



rate proportional to Z2 = Zi + 1 during a time needed to lose the second electron, Tph,2- Applying 



Eq. (|3JJ) repeatedly for all subsequent ionization steps, we get the momentum when the ionization 
is complete. If the total ionization time, is smaller than the Sedov time, in order to get the 
maximum momentum we need to add the further acceleration during the time {tsT — '^'on)' ^i*^ ^ 
rate oc Za?. The final expression for the maximum momentum can be written as follows: 



Pm] 



+ 



^ ZkBkU^TphJc „ , tot \ ^ 

...p. 1.7 z^ +e(^^^-ULT:7^' ^^-^^ 

The last term has been written as a sum over M time-steps in order to handle the case where 
magnetic field and shock speed change with time. 

Let us analyze first the stationary case, with Ush and B constant during all the free expansion 
phase. We consider two different situations which can represent a type I/a and a core-collapse 
supernova (CC SNa). Type I/a SNe explode in the regular interstellar- medium whose typical density 
and temperature are = 1 cm^^ and To = lO'* K. Conversely SNRs generated by CC SNe expand 
into a diluted and hot bubble generated either by the progenitor's wind or by the explosion of other 
close SNe (so called super-bubble)^^. In both cases typical values for the density and temperature 
inside the bubble are = 10^ cm and To = 10^ K. For both type I/a and CC SNe we assume the 
same value for the explosion energy Esn = 10^^ erg, and mass ejecta Mej = \AMq. The resulting 
Sedov times are t^T = 470 yr and t^T =2185 yr, respectively. The average shock speed during the 
free expansion phase is Ush — 6000 km/s in both cases. 

The last parameter we need to estimate is the magnetic field strength which can be inferred 
assuming that both kind of SNe are able to accelerate protons up to the knee energy, which is 
£^knee = 3 • 10^^ cV. This Condition gives By = IbQpLG and 35/iG for type I/a and CC cases, respec- 
tively. It is worth noting that the chosen values of B\ are consistent with those predicted by the 
CR-induced magnetic filed amplification. 

The effectiveness of ionization dependents on the ISRF which is a decreasing function of the 
distance from the Galactic Center. We adopt the ISRF as computed in [14|. In Fig. ([I]) we plot the 



maximum energy achieved for t = tsT by different chemical species, from H up to Zn (Zzn = 30). 
The panels (a) and (b) show the case of type I/a and CC SNRs, respectively. Each panel contains 
four lines: thin solid lines are the maximum energy achieved by ions which start the acceleration 



completely stripped, E^^^, while the remaining lines show ^max computed according to Eq. (3.2) 
for three different locations of the SNR in the Galactic plane: in the Galaxy Center and at 4 and 12 
kpc away from it. Looking at the panel (a) we see that the maximum energy achieved by different 
nuclei in type I/a SNRs does not increase linearly with Z^, instead it reaches a plateau for Zj^ > 25 
at a distance d = 4 kpc and for Z^r > 15 at J = 12 kpc. Only for SNRs located in the Galactic bulge 
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(a) Type I/a: stationary (b) Core-collapse: stationary 




(c) Type I/a: evolution (d) Core-collapse: evolution 



Figure 1: Maximum energy achieved by chemical specie with nuclear charge Zn (ranging from H to Zn) 
at the beginning of the Sedov phase. Panels (a) and (b) show the cases of type I/a and core-collapse SNR 
for a constant shock speed and magnetic field strength, while for panels (c) and (d) Msh and B evolve with 
time, as explained in the text. The thin solid line shows the maximum energy achieved by atoms which 
are fully ionized since the beginning of the acceleration, while the other curves are computed including the 
photo-ionization due to ISRF for SNRs located at three different distance from the Galactic Center, namely 
0, 4 an 12 kpc. 



the proportionality relation ^max °^ '^n holds, at least for elements up to = 30. The effect of 
ionization is much less relevant for core collapse SNRs: only those remnants located at a distance 
of 12 kpc show a noticeable reduction of ^max- 

Previous conclusions are based on the assumption that both the magnetic field and the shock 
velocity remain constant. This is indeed a poor approximation. In fact during the free expansion 
phase the shock speed can vary by a factor of few. Now for DSA the acceleration rate is propor- 
tional to uli^5B{ush), where 5B{ush) is the turbulent magnetic field generated by the CR-induced 
instabilities (resonant or non-resonant), which is, in turn, an increasing function of m.,/,. Hence the 
acceleration rate can significantly change during the free expansion phase. In order to evaluate the 
effect of the time evolution we can use a simple toy model always in the framework of steady-state 
approach: we approximate the continuum evolution into time steps equal to the ionization times, 
Tph,k> assuming that the stationary approximation is valid during each time step. For the description 
of the shock dynamic we follow [jT5|]. Specifically we adopt the solution for a remnant characterized 
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by a power-law profile of the ejecta with index n = l, expanding into an homogeneous medium. For 
this specific case the shock velocity for t < tsr is well described by u\{t)/uch = 0.606 (f/fc/i)^^^^ 
(see Table 7 in ^) where Uch = {Esn /Mej^l'^ , tch = E^^'^mI^ p^^'^ and tsT = 0.132t,.h. We 
adopt the same values of Esn, Mej and pi used in the cases (a) and (b). 

For what concern the magnetic field we assume that the amplification mechanism converts 
a fraction of the incoming kinetic energy flux into magnetic energy density downstream of the 
shock: 5B^{t)/{S7i) = aBPiuf{t). The parameter hides all the complex physics of magnetic 
amplification and particle acceleration. For the sake of completeness we mention that in the case 
of resonant streaming instability as °^ ^ckVa/ui, where ^^r is the efficiency in CRs and va is the 
Alfven velocity, while in the case of resonant amplification oc ^^i.ui/4c [^. However such 
relations cannot be applied in a straightforward way since they require using a non linear theory. 
Here we prefer to make the simplest assumption, taking ag as a constant. As we done for the 
stationary case, we determine the value of as assuming that for t = tsr the energy achieved by 
protons is 3 • lO'^ eV. This condition gives ttg = 3.25 • 10"^ and 6.80 • 10"^ for type I/a and CC 
cases, respectively. It is worth noting that in the case of some young SNRs the value of as has 
been estimated from the measurement of both the shock speed and the magnetic field strength and 
the results are only a factor 5 — 10 larger then the values we use here (see, e.g.. Table 1 from [^). 

Now Eq. (3^) can be used to evaluate the effect of time evolution. For each time-step, Tph./t, 
the values of and are computed according to the equation described above, evaluated at the 
beginning of the time-step. Always in Fig. |I| we report the results for the maximum energy att = tsr 
for the case of type I/a and CC SNRs in panels (c) and (d), respectively. Our results show that 
the reduction of the maximum energy is now more pronounced with respect to the stationary case 
shown in the panels (a) and (b). Even in the scenario of a core-collapse SNR located in the Galactic 
bulge, iron nuclei are accelerate only up to about one half of the maximum theoretical energy. The 
reason why ionization is less effective when the evolution is taken into account is a consequence 
of the fact that acceleration occurs mainly during the first stage of the SNR expansion. In fact both 
the shock speed and the magnetic field strength are larger at smaller times. On the other hand the 
effective ions' charge is small during the initial phase of the expansion, hence the acceleration rate 
is smaller than its maximum possible value. 



4. Discussion and conclusions 



The most relevant consequence of the ionization mechanism concerns the shape of the knee 
in the CR spectrum. As we have already discussed, the relation f'max.A' °^ is needed in order to 
fit the data in the knee region, which has a slope equal to 3.1. Even a small deviation of the cutoff 
energy from the direct proportionality can affect this prediction. Our results suggest that type I/a 
SNRs seem unable to accelerate ions up to an energy Z^r times the proton energy, due to their small 
Sedov-Taylor age. This consideration is strengthened by the fact that the CR flux observed at the 
Earth is mostly due to the SNRs located in the solar neighborhood, rather than those located in the 
Galactic bulge' where the larger photon field make the ionization faster and allow nuclei to achieve 

' In fact the escaping length from the Galaxy is determined by the thickness of the Galactic halo which is --^ 3 — 5 
kpc, hence particles can reach the Earth only if they are produced within this distance, while the Galactic bulge is at 8 
kpc from us. 
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larger energies. Ions accelerated at SNRs produced by CC SNe are less affected by the ionization 
and can achieve larger energies because the larger Sedov-Taylor age. But even in this case the 
maximum energy could be appreciably reduced when the time evolution of the remnant is taken 
into account. We can conclude that the primary sources of CRs above the knee energy are most 
probably the core-collapse SNRs with Mej » \Mq. 

A second comment concerns the acceleration of elements beyond the iron group. Even if the 
contribution of such ultra-heavy elements to the CR spectrum is totally negligible at low energies, 
in the 100 PeV regime it could be significant. This region is where the transition between Galactic 
and extragalactic CRs occurs. Indeed several authors pointed out that a new component is needed 
to fit this transition region, beyond the elements up to iron accelerated in "standard" SNRs (see 



e.g. the discussion in [Q]). As inferred by [ [rT| ] stable elements heavier than iron can significantly 
contribute to the CR spectrum in the 100 PeV regime if one assume that their maximum energy 
scales like Z^. This assumption is especially appealing also because it could explain the presence 
of the second knee in the CR spectrum [[IT]]. In fact the ratio £'2<.<;^„gg/£'fa,ee ^ 90 is very close to the 
nuclear charge of uranium, which should have ^max.f/ = 92£'niax,p = 414 PeV. On the other hand 
the ionization mechanism discussed here provides a strong constraint on the role of ultra-heavy 
elements. If the acceleration of elements heavier than Fe occurs at SNRs like those considered 
in this work, it is easy to show that they cannot achieve energies much larger than the Fe itself, 
because the total ionization time increases rapidly with the nuclear charge. Even if we consider 
the acceleration in very massive SNRs, with Mgj » \Mq, the contribution of ultra-heavy elements 
remains unlikely. An interesting alternative worth to investigate is the acceleration during the very 
initial stage of the explosion, when the expansion occurs into the dense progenitor's wind and the 
Coulomb collisions can strongly reduce the ionization time. 
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